Crossover behavior for complex order parameter in high- Tc superconductors 
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A number of recent experiments have suggested the presence of either real or complex compo- 
nents in the gap symmetry of high- Tc superconductors (HTSC). In this paper we introduce a novel 
approach to study the competition of such complex order parameter mixtures by varying the position 
of the two-body attractive potential in a two dimensional extended Hubbard Hamiltonian. We show 
that this procedure explain a number of experimental results and on the theoretical side, it may be 
related with certain HTSC microscopic models like the spin fluctuation theory. Following current 
trends we concentrate on the study of d^2_y2 order parameter with a component of the type d^y 
or a s-wave like s^2^y2 and Sxy symmetry. We show that the position of the optimal s-component 
peak changes with the position parameter b while the d-component occurs always in the optimally 
region around hole content p ^ 0.39. These studies may be useful to interpret some experimental 
data and to explain why similar experiments yield different gap symmetries. 
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I. INTRODUCTION 



In order to understand the fundamental mechanism 
for high- Tc superconductors (HTSC), many experiments 
have attempted to find out the properties of the super- 
conductors pair wave function. Notwithstanding this ef- 
fort the nature of the orbital order parameter symmetry 
(or energy gap) has not yet been settled, despite increas- 
ing evidences toward a major d^^-y^ symmetry as it ex- 
plains a number of different experimental results 
Furthermore the d-wave state has also theoretical sup- 
port from weak and strong-coupling approaches |^] and 
from Monte Carlo and numerical studies on a two di- 
mensional Hubbard model Q. On the other hand sev- 
eral experiments have also suggested the presence of real 
or complex mixture of order parameters. For instance, 
c-axis Josephson tunneling data between twinned YBCO 
and a s-wave superconductor were interpreted by a con- 
densate containing a mixture of dx^-y^ and s order pa- 
rameter 1^,^. This type of order parameter symmetry 
was argued to be a consequence of the orthorhombicity 
of the YBCO crystal structure [^|j8j . Recently, Klemm at 
al 1^] have studied the symmetry operation of the crystal 
groups relevant to HglOOl, YBCO and BSCCO to ana- 
lyze the possibility of different order parameter symme- 
try combinations for these HTSC and their results agree 
with the c-axis Josephson tunneling experiments. Fur- 
thermore the possibility of complex mixtures of differ- 
ent symmetries has also been suggested as it would 
induce time reversal symmetry breaking in connection 
with magnetic defects or small fractions of a flux quan- 
tum <i>o = hc/2e observed in YBCO powders. Similarly 
a complex admixture of (i^2_j,2 and s symmetry was also 
proposed to explain the ab-oriented YBCO/I/Cu 
tunnel junctions data ]T2[ |. They proposed an interest- 
ing and very simple physical picture: Andreev scatter- 



ing near the YBCO surface causes strong pair dissoci- 
ation. The resulting quasiparticles may then be paired 
again also by a subdominant s-channcl that is less sensi- 
tive to surface pair breaking then the dominant dx^^y^- 
channel. Calculations minimizing the free energy show 
JlH that such complex mixture can coexist at low tem- 
peratures. Using a tight binding model in a mean field 
treatment Honerkamp et al |13| have arrived at a simi- 
lar conclusion that a surface phase transition towards a 
time-reversal symmetry breaking surface state carrying 
an s+idx2^y2 state may appear (below Tc). Recently 
a complex mixture of d^2_j,2 -I- idxy was also proposed in 
order to explain a phase transition characterized by mag- 
netic field plateaus or kinks in the thermal conductivity 
|]l5| of Bi2Sr2CaCu20s compounds. 

Based on the above discussion, we study in this paper 
the competition of the dx2-y2 and others subdominants 
complex symmetries of the gap function. We develop in 
these calculations a novel approach based on a change 
in the position of the attractive potential V from its 
usual nearest neighbor position in the extended Hubbard 
Hamiltonian. As it is well known, the on-site Coulomb 
correlations may explain the antiferromagnetism of the 
low doping regime, the large magnetic fluctuations and 
the semiconductor like properties of the metallic phase 
|]l6|-^^. Furthermore the presence of a small attractive 
interaction leads to phase separation and, on a square lat- 
tice, yields a superconducting phase with a dx2-y2 sym- 
metry 1^. Thus, the Hubbard Hamiltonian with such 
phenomcnological attraction is a natural candidate to 
deal with the HTSC. The idea to change the position of 
the attractive potential in the extended Hubbard Hamil- 
tonian was originally motivated by the experimental fact 
that different compounds have different Tc x n phase di- 
agram and the optimal doping values may vary accord- 
ing the compound. This fact can be physically inter- 
preted due to a possible change in the range of the attrac- 
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tive interaction that leads to pair formation for 
compounds with short range correlations the carrier den- 
sity must be larger than those compounds with a longer 
ranged interaction in order to undergo a superconduct- 
ing phase transition. The same type of physical argument 
can be used to interpret the measured different values of 
the coherence length ^ for different family of compounds 
since a short range interaction requires larger densities 
than a long range one in order to produce the coherence 
motion that leads to superconductivity. Furthermore, 
the typical parameter which characterizes the strength 
of the attractive interaction |2^, x = rneO^V /2fi'^ ^ when 
estimated for most of HTSC, varies between 10-3000 and 
such high values indicate that the size of the bound states 
(or Cooper pairs) are indeed related to the size of the 
minimum of the attractive potential. 

Another interesting aspect from the theoretical point 
of view is that this parametric change of the attractive 
potential can be related with the order parameter ex- 
pansion introduced by the spin- fluctuation theory [p2| . 
In this approach it has been proposed a d2;2_j^2-wave gap 
of the form AQ{k)(coskxa — coskyo) where Ao(fc) is ex- 
panded in powers of {coskxO, -\- coskyO). Taking a close 
look at this expansion, one can easily verify that it con- 
tains terms like {coskxa—coskya){coskxa+coskya) which 
are proportional to {cos(2kxa) — cos{2kya)) and what can 
be seen as a type of d-wave gap symmetry that arises 
from a potential like V{k) — VQ{cos{2kxa) + cos{2kya)). 
By the same token, we can find terms proportional 
to {cos{3kxa) — cos{ikya)) which can be originated by 
a potential V{k) = Vo{cos{3kxa) + cos{Skya)) and so 
on. There will also be crossed terms proportional to 
(coskxCL X coskyo) that may be associated to a next- 
nearest potential. Thus, we see that the gap expres- 
sion proposed by the spin-fluctuation theory contains 
several terms and some may be associated to a poten- 
tial of the form V{k) — Vo{cos{bkxa) + cos{bkya)), where 
b = 1,2,3,... is a parameter related to the position of 
the attractive potential (6=1 leads to the usual nearest- 
neighbor expression). Consequently we have studied in 
this paper the phase diagram of complex admixtures of 
a d-wave like order parameter with some minor complex 
component with this parameterized change of the poten- 
tial position, that is, with different values of b. These 
calculations are made with the extended Hubbard model 
on a square lattice of side a, which describes a tight bind- 
ing model with the strong correlations taken into account 
and is defined by 

H = - ^iii^tcja + h.c.) + UY^ n.iTii^ - V^niTij. 

{i]),<y i (ij) 

(1) 

where is the hopping integral between sites i and j 
and are estimated by comparison with either band struc- 
ture calculations or Fermi surface measurements, U is the 



on-site correlated repulsion and V is an attractive phe- 
nomenological interaction which in principle can be due 
to spin fluctuations, as mentioned above or any other 
type of mechanism. In order to obtain quantitative re- 
sults, we have used the tight-binding dispersion relation 
derived from Fermi surface measurements by Norman et 
al for Bi2Sr2CaCu20s but this choice is not a funda- 
mental requirement, another studies have used different 
dispersion with small quantitative changes. 

II. THE METHOD 

We used a BSC mean-field approximation ||2^,|2^ to 
study the superconducting phase associated with the 
Hamiltonian of Eq. (|l|) described above. In this model 
the superconducting state is characterized by a gap order 
parameter, which at a finite temperature T satisfies the 
BCS self-consistent equation 

A^--E%^r (2) 
r 

where Vjj- is the Fourier transform of the potential of Eq. 
® and 

= ^*""^2St ' ^ V^(^r^^^^^)M^ (3) 

where is the dispersion relation taken from Ref. 
p3[ and /j, is the chemical potential. Following along 
the lines of Ref. , one can easily show that Vj^j- 
may be written in a "separable" form, Vj:f ~ U — 
2V {cos{kxa) cos{lxa) + cos(fcya) cos{lya)) which leads to 
the usual d and s-wave. This procedure can easily be gen- 
eralized for the case of 6 = 2, 3... discussed above and will 
be developed below in order to study the superconducting 
transition with respect to the position of the attractive 
potential. For a two-component order parameter with a 
separable form, the corresponding gap equations may be 
written as 

%-E^i1^'% = EM (4) 

Here, since we want to study the competition be- 
tween the dx2-y2 and others symmetries, we have 
always taken Vi — Vd^2_y2i f\ = (cos(bkxa) — 
cos(bkya)) and Ai ~ Ad ^ 2 and the complex com- 
ponent can either be V2 = ^xyi^s^2+y2 o^' ^Sxy with 
/| = 2sin{bkxa)sin{bkya), {cos{bkxa) -\- coslbkya)) or 
2cos{bkxa)cos(bkya) and the gap functions are Ai = 
iAd^y, As^2+ 2 or A^^j^ respectively. Now, separating the 
real and imaginary parts and combining Eq.|| and Eq.^ 
we obtain the following two gap equations 
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These two equations has to be solved with the density of 
carriers p equation which is given by [Eq] 
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Thus we solve self-consistently the above 3 equations 
(Eq.|| and Eq.||) in order to study the phase diagram 
of a complex admixture of the order parameter. These 
equations are general in the sense that they can be ap- 
plied to any type of two component symmetry mixtures. 
The results will be discussed below. 



III. RESULTS AND DISCUSSIONS 
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FIG. 1. The critical temperatures for the s and d-channels 
for 6=1. The continuous line is the calculations described 
in the paper and the diamonds are the experimental points 
from Ref. 28 

We start with the analysis for the complex mixing 
with the two most relevant ^ symmetries, The dj;2_y2 
and the Sx2^y2 which we write as d^^-y^ + is^^+yi. We 
have solved the above self consistent equations with fixed 
Vd/t — 0.11 and Vs 2_ 2/* = 0.18, which these values 
were chosen in order to obtain reasonable values for the 
critical temperatures Tc in the range of parameters stud- 
ied. In this way we obtain a maximum Tc ~ 94K for b=l 
and by the dominant d^2_y2 symmetry near the optimum 
doping. In Fig.l we show the results for the values of 
Tc X p in order to compare with the experimental results 
of Allgeier et al |2|] on Bi2Sr2CaCu20s+5 (Bi2212). No- 
tice that their maximum Tc occur at p w 0.32 while our 
calculations yield a maximum Tc at p ~ 0.39 and this 



discrepancy is probably due the dispersion relation that 
we used since a similar BCS mean field calculation using 
the Hubbard Hamiltonian (for b=l) gives a maximum 
which agrees with the experimental result |24|. 



We should point out that the comparison with a HTSC 
experimental phase diagram, as in Fig.l, has only mean- 
ing if the material has a BCS type behavior, that is, when 
the appearance of the superconducting gap and phase 
coherence occur simultaneously at the onset of super- 
conductivity. On the other hand, angle-resolved photoe- 
mission spectroscopy (ARPES) on underdoped Bi2212 
has revealed the presence of a "pseudogap" above Tc 
p9| , ^ 23 1 which has also been confirmed by recently 
Fermi surface ARPES measurements This unusual 
behavior may be an indication of non-BCS behavior 
which Tc is controlled by the doping level, at underdoped 
regime, rather then by phase coherence with a single par- 
ticle gap ||3^. Nevertheless this is not a completely 
settled scenario which only more experiments will give 
the correct interpretation specially because the pseudo- 
gap has not yet been seen in several cuprate supercon- 
ductors Isl. 
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FIG. 2. Results for the zero temperature A^(^^2_y2-^ (con- 
tinuous lines) and Ag(^^2^y2^ (dashed lines) in units of 
t = 0.26eV as function of the hole density p for 6=l(a),2(b) 
and 3(a). p = is the half-filling density. 

In the Fig. 2 we show the results for the zero tem- 
perature gap amplitudes ^ ^ (continuous line) and 
As .-.^ 2 (dashed line) as function of the density p for 
three values of the parameter b, namely b — 1,2 and 3. 
Where there is superconductivity region, the values of A 
have their maxima at zero temperature and vanishes at 
Tc in agreement with the BCS method. In fact, the self- 
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consistent calculations yields the temperature, the den- 
sity p and the value of the gap as one solves Eqs. 5 and 
6 simultaneously. In Fig. 3 we show how the gap changes 
with the temperatue. The zero temperature gap value 
is proportional to Tc and the constant of proportionality 
differs slightly from one channel to the other. We get 
/ksTc « 2.6 in the optimally region. 
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FIG. 3. The display of the doubled crossover. As the 
density increases we see the solutions change from dominant 
s-channel in a), to a dominant d-channel as in b) and again a 
dominant s-channel in c). The temperature T is in Kelvins. 



There are very interesting differences from the s and 
d channels which arises as we vary the parameter h. For 
6 = 2, the values of Ad^2_ 2 are basically suppressed 
while for 5 = 3 its its strength is decreased by about 
half of its value for 6=1 (which agrees with the calcu- 
lations of Ref . ) and does not move appreciably from 
the optimum doping value of hole content p — 0.39 (p is 
the fraction of holes per Cu atom in the CUO2 sheet). 
This result is a little higher then the experimental result 
p8| which yields p = 0.32 for the Bi optimum doping 
value. The maximum Tc is related with the envelope of 
the curve shown in Fig. 2. On the other hand, for the 
extended s, Ac , , , the change in 6 causes the value 
of the corresponding optimum doping density to change 
continuously; it is at p « 0.41 for 6 = 2 and and about 
p = 0.32 for 6 = 3 as it is seen in Figs. 2a, 2b and 2c. At 
6=1 (Fig. 2a), the s-wave has the optimum doping at 
p ~ 0.05, for 6 = 2 the maximum occurs at p « 0.40 
and it dominates over A^; which is almost entirely sup- 
pressed. For 6 = 3 the competition between the two 
channels becomes very interesting; the maximum of As 
has moved to p Ri 0.33 while the maximum of A^ re- 



mains fixed at p « 0.42 which gives rise to a completely 
novel doubled crossover phenomena; up to p « 0.35 the 
self-consistent equations yields pure extended s-wave so- 
lutions. Between p « 0.35 — 0.37 the mixing of order 
parameters starts at low temperatures but with the s- 
channel dominating over the d-channel as its critical tem- 
perature Tc is larger as shown in Fig. 2c and at Fig. 3a for 
p = 0.37. 

At p « 0.38 a crossover arises and the d-channel so- 
lutions dominates over the s-channel which becomes the 
minor component as shown in Fig. 2c and Fig. 3b. There is 
always mixing at low temperatures (the continuous curve 
of the minor s-wave never vanishes completely) but the 
d-channel is the only one at Tc- At p w 0.48 there is a 
new crossover and the dominant solution changes again 
as the s-wave becomes the dominant solution and the 
only one at Tc and the d-wave appears again at low tem- 
peratures as the minor components as displayed in Fig. 2c 
and Fig. 3c. Notice that the appearance of this subdom- 
inant component is exactly as the mechanism proposed 
to explain the ab-oriented YBCO/I/Cu tunnel junction 
data 1 12 1 discussed in the introduction. 
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FIG. 4. Calculations for the A^(^^2_y2^ (continuous lines) 
and As(2,y) (dashed lines) in units oi t = 0.26eV as function 
of the hole density p for 6=l(a),2(b) and 3(a). 



Next we have studied the complex gap function of sym- 
metry d^2_y2+isxy The results are quite similar to those 
above; for 6=1 the d-state dominates and has its opti- 
mum doping at p = 0.40 as shown in Fig. 4a. For 6 = 2 
the d-channel is again suppressed and the Sa;j,-state is the 
only one and its maximum doping is also at p ~ 0.83 as 
displayed in Fig. 3b. For 6 = 3 there is again the very 
interesting doubled crossover between these two compet- 
ing symmetries. For 0.2 < p < 0.36 there is only s-wave 
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instability but at p « 0.36 a small d-wave channel devel- 
ops at low temperatures. Between p k, 0.38 and p ~ 0.45 
we find only pure (1^2 _y2 solutions as opposed to the case 
shown above in Fig. 2c where a low temperature mixing 
was always present. For p > 0.45 the s^^y-channel starts 
to develop, initially as a minor component (see Fig. 5b) 
and above p — 0.47 a new crossover appears as it be- 
comes the dominant channel and the only channel at Tc- 
Thus we verified that the states dx^-y-i + iSxU has also 
the doubled crossover and the dominant regions are sim- 
ilar to those for dx2-y2 + isx2^y2 however, there is a new 
feature with this type of admixture: we can see from 
Fig. 4c that, for 6 = 3, one kind of gap symmetry influ- 
ences the other. One can easily see the abrupt change 
in As^^ near p w 0.33 and the change in slope for both 
channels near p « 0.48. The new phenomenon, that is, 
the slightly suppression on A(T) of the dominant chan- 
nel can also be seen on Figs. 5a, 5b and 5c. This is a 
new feature not seen in Figs. 3c, and this suppression is 
the cause of the pure d-wave solutions at the optimal re- 
gion which did not occur for previous studied extended 
s-wave. This interesting feature may be similar to what 
happens in Andreev scattering near surfaces as discussed 
in the introduction. 
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around the same optimal region and the other channel 
gives only a very small component at low temperatures. 
Thus we do not observe any crossover for dj.2_y2 -|- id^y 
complex mixture, simply the interchange of dominant 
symmetry as we change b. 



IV. CONCLUSIONS 

The complex mixtures in superconducting order pa- 
rameter were studied here as function of hole density in 
the context of the BCS mean-field calculations for the 
extended Hubbard model. We have revealed very in- 
teresting properties that may be relevant for HTSC; if 
the interaction is sufficiently short ranged that can be 
mapped by a nearest neighbor potential, at low doping, 
the system is described by either pure dj.2-y2 or Sj;2+y2 
order parameter. If the range of the potential is larger the 
system changes its order parameter symmetry (for 6 = 2) 
or may admit a mixed state with both s and d compo- 
nents present (for 6 = 3). In this last case the dominant 
component depends strongly on the density or doping 
level p and may change abruptly the dominant symme- 
try at Tc- Since in most HTSC it is difficult to determine 
p accurately ]35| and as our calculations demonstrate, a 
small change in p may be accompanied by a crossover 
of the dominant component. If we assume an expansion 
with 6 = 1,2,3... like that of the spin fluctuation theory 
of Ref. ||2^, then one must look at their weight to see 
which one is more important and concomitantly, which 
symmetry of the order parameter is the dominant one. 
This present work provides a hint to understand why 
different order parameter symmetries appear in different 
experiments performed on the same type of compounds. 
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FIG. 5. The display of the doubled crossover. As the 
density increases we see the solutions change from dominant 
s-channel in a) , to a dominant d-channel as in b) and again a 
dominant s-channel in c). The temperature T is in Kelvins. 



Lastly we have also studied the complex order parame- 
ter mixing dx2-y2 +idxy for the three values of b, however 
we have not found even a single crossover. For 6=1, 
the o?2,2_y2 channel dominates and has also its optimum 
value at p = 0.40 and the dxy appears only at low tem- 
peratures. For 6 = 2 and 3 the dxy dominates completely 
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